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To probe the nature of the even-fraction quantum Hall system, we have investigated the low- 
lying excitation spectrum by means of exact diagonalization for finite systems. We have found (i) 
a striking one-to-one correspondence (i.e., a shell structure) between the spectrum and those for 
free (composite) fermions, (ii) a surprisingly straight scaling plot for the excitation energy that 
gives a zero gap (metal) in the thermodynamic limit, (iii) the effective mass evaluated from the 
scaling becoming heavier for v = 1/2, 1/4, 1/6, but (iv) some deviations from the single-mode or the 
Hartree-Fock composite fermion approximation as well. 
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In the physics of the fractional quantum Hall system, 
the composite fermion (CF) picture H not only serves 
as an illuminating way of understanding Laughlin's in- 
compressible quantum liquid for the odd-fraction Landau 
level filling, v, but also poses an interesting question of 
what is the nature at even fractions, which is the accu- 
mulation point of the fractional quantization. A seminal 
paper by Halperin, Lee, and Read Q suggested that the 
system at v = 1/2 should be a Fermi liquid of CF's in 
the mean-field picture, which led to intensive studies. In 
contrast to the incompressible quantum Hall state or su- 
perconductors where energy gaps arise from many-body 
effects, we have to question here how the gap vanishes 
(i.e., how the liquid becomes compressible) despite the 
presence of the electron correlation. 

Naively a CF, composed of an electron and an even 
number (<f> = 2, 4, • • ■) of flux quanta, feels the mean mag- 
netic field B c g = — 4>)4>op, where B = v~ l (j)Qp is the 
external magnetic field, p the number density of elec- 
trons, and 0o = 27r/e the flux quantum (in the units 
in which c = 1 and h = 1). Thus B c g vanishes for 
v = l/(f>. There is, however, no guarantee that the mean 
field should be good, and the above argument does not 
in fact say anything as to where the electron-electron 
repulsion comes in. Recent developments [^[jH) have in- 
dicated that we can define a 'dipole' (composite parti- 
cle + a correlation hole), where the flux-attachment is 
thought to mimic the repulsive correlation of electrons. 
The Halperin-Lee-Read prediction on v = 1/2 has been 
re-examined in the dipole picture, and the compressible 
nature is reproduced M. 

These approaches still adopt mean-field treatments, 
and their validity has yet to be fully clarified. The dif- 
ficulty arises because fluctuations of the Chern-Simons 
gauge field that implements the flux-attachment should 
be significant. The fluctuations in fact determine the 



residual interaction between CF's as well as the effective 
mass, m*, of a CF, which are difficult to evaluate an- 
alytically. Hence exact numerical studies for finite sys- 
tems are valuable. Rezayi and Read have numerically 
shown that the ground state for the v = 1/2 system on a 
sphere has the same angular momentum as expected from 
Hund's second rule for the same number of fermions in 
B = 0. Morf and d'Ambrumenil Q have estimated m* 
from the size scaling of the ground-state energy. How- 
ever, we are still some way from understanding to what 
extent the CF picture applies. 

One direct way going beyond the ground state is to 
look at the low-lying excitation spectrum — here we 
question whether or not there is a one-to-one correspon- 
dence, in the structure of the excitation spectrum, be- 
tween the v = 1/2 liquid and a free fermion system in 
B = 0. This can also enable us to extract, thorough the 
size-scaling of the energy gap, the effective mass. This is 
exactly the motivation of the present work. 

There are two points we wish to make: (i) how to 
perform the size scaling is always a subtle problem, es- 
pecially so when detecting the excitation gap that may 
vanish in the thermodynamic limit, (ii) some analytic 
studies HH have indicated that the nature of the v = 1/2 
liquid is affected by the range of the electron-electron in- 
teraction. So we have taken a specific scaling sequence, 
and also varied the range in monitoring the excitation 
spectrum. 

We shall show that, (i) we do have a striking one-to-one 
correspondence between the interacting and free systems. 
The shell structure in the spectrum is deformed with 
the range of the interaction, which is interpreted here in 
terms of the single-mode approximation (SMA), (ii) the 
effective mass becomes heavier as v = 1/2 — > 1/4— > 1/6, 
somewhat more slowly than the Hartree-Fock (HF) pre- 
diction of mjfjp oc 1/v 2 . 
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We adopt the edge-free spherical system following Hal- 
dane, which has an extra virtue that the full rotational 
symmetry can be exploited in classifying the states. 
Dirac's quantization condition dictates that the total flux 
AirR 2 B be an integral (25) multiple of the flux quan- 
tum, where R is the radius of the sphere. The eigen- 
value of the non-interacting part of the Hamiltonian is 
e = [1(1 + 1) - S 2 }/(2mR 2 ), where I (> S) is an integer. 
The lowest Landau level (LLL) corresponds to I = 5 with 
the Landau level filling given by v — (N — l)/25 for N 
electrons. 

For later reference let us look at the spectrum for free 
fermions in zero magnetic field, which is the mean-field 
solution for CF's at v = 1/2. The energy of a free fermion 
on a sphere is readily given by e — 1(1 + l)/(2m*R 2 ), 
where l(> 0) is the angular momentum and m* the 
fermion's mass. We have to note that, with each level 
being (21 + l)-fold degenerate, a 'closed shell' configura- 
tion is realized when N — (lp + l) 2 [Fig. 0(a)]. Here 
If(= 1, 2, 3, • • • for N = 4, 9, 16, • • •) is the highest occu- 
pied I, so that we may call this the Fermi angular momen- 
tum in analogy with the Fermi momentum in the planar 
system. 
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FIG. 1. (a) A closed-shell ground state of the N = 16 
free system. Solid (open) circles represent occupied 
(empty) states, (b) An example of multi-exciton excitations 
([l F -l][l F ] 2 -» [Zf + 1] 3 here). 

When TV ^ (lp + l) 2 , the ground state of the non- 
interacting system is thus degenerate, or has an 'open 
shell'. For interacting particles the degeneracy is lifted, 
and the total angular momentum of the ground state be- 
comes nonzero 0. Since this can obscure the scaling, a 
more straightforward way is to concentrate on the closed- 
shell sequence satisfying N = (lp + 1) 2 =4, 9, 16, • • •. For 
this sequence the total angular momentum of the ground 
state remains zero, and provides a natural sequence to- 
ward the infinite system for establishing both the struc- 
ture of the low-lying excitation spectrum and the energy 
gap. 

The simplest class of excitations from a closed shell 
is 'single-exciton' excitations where a particle is ejected 
from the Ip-th shell to the (lp + l)-th, as has been pointed 
out by Rezayi and Read Q|. The exciton's angular mo- 
mentum takes the values L = 1, • • • , 21 p + 1. These exci- 
tations (abbreviated here as [lp] — ► [lp + 1]) provide the 
lowest-lying branch for 1 < L < 2lp + 1. 

We can generalize this, including multiple excitons, to 
obtain the whole picture. For 2lp + 1 < L < Alp, the 



lowest-lying excitations are [If] 2 — > [If + I] 2 , i-e., two- 
exciton excitations. For Alp < L < 61f ~ 3 for N > 9 
[lp] 3 — > [lp + l] 3 and so on, where n-exciton excita- 
tions [lp] n -> [I p + 1]" exist for L < n (2l F + 2 - n). 
For L > 61 p — 2, more complicated excitations such as 
[; F -l] m [; F ] n -> [Z F + l] m '[; F +2] n ' (m+n = m'+n') must 
also be considered. Overall, however, the lowest-lying 
states are one-, two-, three-, • • • excitons, whose energies 
A e form steps moving up at Lmax — 21f, Alp, 6Z F , • • •, re- 
spectively, as shown in Fig. ||for N = 16, although there 
are finite-size corrections in -Lmax = 2lp + l,Alp,Qlp — 3, 
etc as we have seen. 

Having looked at the free case, we now come to the 
structure of low-lying excitations in the interacting sys- 
tem. The exact low-lying energies are obtained by diag- 
onalizing the Hamiltonian matrix. For v = 1/2 we have 
25 = 2(./V — 1), and the dimension of the Hamiltonian is 
4,669,367 in the L z = subspace for N = 16 electrons. 
The matrix elements can be expressed in terms of Hal- 
dane's LLL projected pseudo-potential, Vj Jl(],[ll]]. If we 
explore the evolution of the spectra to N = 16 (Fig. 0) 
from those for N = A (not shown) and N — 9 (Fig.^5| 
below), we are led to a well-defined series of cusps in the 
excitation spectrum, whose positions exactly agree with 
the above-mentioned positions for the free fermions. The 
degeneracies in the latter case are naturally lifted due to 
the interaction between CF's, but the interaction is weak 
enough to preserve the shell structure, which remarkably 
persists up to the angular momentum as large as 30 [jl3| . 
This is the first key finding in this Letter. 




FIG. 2. Low-lying excitation spectrum for N — 16 (solid 
circles). Dashed line is a guide for the eye, while the full curve 
represents L(L + 1). The low-lying excitation spectrum A e 
for free fermions is also shown (□). 

We can next evaluate the energy gap, A. In the 
free fermion system, the lowest excitation corresponds 
to lp ->• lp + 1 with A ee (l F + l)/(m*R 2 ). This quantity 
has a well-defined scaling, A = (An p / m*)[\/N / (N - 1)] 
when N is varied with p fixed. For the interacting system, 
the cusped structure revealed here enables us to identify 
the position of the lowest excited state, which always oc- 
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curs for the first cusp at L = 21 p + 1 (the high-L end of 
the single-exciton excitation). Fig. || shows this gap for 
v = 1/2 We can immediately see a surprisingly ac- 
curate linear scaling that extrapolates to zero for N — > oo 
if we have y/~N /(N — 1) as the scaling variable, as guided 
by the free-system behavior, A = (4ivp/m*)[\/~N / (N—l)] 

N 
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FIG. 3. Size scaling of the gap for v = 1/2, 1/4, 1/6. The 
dashed lines are linear fit to the data. 

This same formula can be used to extract the effective 
mass m* of CF's from the gradient of the scaling plot, 
with the result 1/m* = (0.185 ± 0.002)e 2 ^, where I = 
1/yeB is the magnetic length. The 1/m* obtained here 
from the excitation gap is slightly smaller than 1/m* ~ 
(0.2 ± 0.02)e 2 £, estimated from the ground-state energy 
per particle ||. On the other hand the present value is 
slightly larger than the analytic estimate, 1/m* ~ e 2 £/6, 
obtained from the interaction energy between an electron 
and a correlation-hole in the first-quantized picture || or 
the self-energy of the CF in the temporal gauge in the 
HF approximation Jl5"| ]. 

The gap and mass, dominated by gauge field fluctua- 
tions, should depend on the number of flux quanta at- 
tached ((f)), so we further obtained the scaling plot for 
the sequence v = l/4> = 1/2,1/4,1/6 in Fig. | @. 
The gap again vanishes for N — > oo, where the effec- 
tive mass systematically becomes heavier in the sequence 
v = 1/2, 1/4, 1/6 as shown in Fig. | for N = 9. In the 
HF approximation for the CF we can show |n| that m* 
should scale as 
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This is a decreasing function of <j> as well, but the present 
numerical result is seen to deviate from the HF result 
(inset of Fig. 0) for larger 0. Q 
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FIG. 4. A (oc 1/m*) for the sequence v = 1/2, 1/4, 1/6 for 
TV = 9. The dashes lines represent the HF result, with the 
inset depicting m*j F /m*. 

Now let us look more closely at the excitation spec- 
tra. Note in passing that the overall shape of the 
spectrum exhibits an oc L 2 asymptote as evident 
from Fig. |^. We can explain this by converting the 
Hamiltonian in the c'c'cc form to c^cc^c. We have 
then, up to a constant, e 2 /(2£\f~S) Y^k=i Pk • Pk 
where pa ■ Pk = 1) M Pkm Pk-m and p K M = 

Y J (-) s+m2 (Sm 1 ;Sm 2 \KM) c^c-m, where cre- 
ates the mth orbit. The transformed coefficient be- 

comes V K = EjloM) 2 ^ 2 J + l){ff*r}^ whcrc 
(ji m i', j2iri2\JM} is the Clebsch-Gordan coefficient and 
{flif } Wigner's 6j symbol. In this representation, while 
Pim is nothing but the (LLL projected) total angu- 
lar momentum operator, so the leading term becomes 
pi ■ pi = [3/S{S + l)(2S + l)]L ■ L, which explains the 
asymptote oc L(L + 1). 

Now we come to what happens when the range of 
interaction is changed. We have calculated the excita- 
tion spectra replacing the pseudo-potential V%s-m with 
(Vas— m) a ■ Since ~Vis-m is the potential between two 
electrons with the relative angular momentum m, a < 1 
(a > 1) corresponds to the interaction longer- (shorter- 
ranged than Coulombic. 

The numerical result in Fig. || |l3| shows that the 
cusped structure in the spectra becomes more pro- 
nounced (i.e., effect of the inter-CF interaction becomes 
enhanced) as the interaction is made shorter-ranged, al- 
though the positions of cusps remain the same. So the 
free CF picture seems to be better for longer-ranged in- 
teraction. This is in sharp contrast with the Laughlin's 
liquid at odd denominators for which the mean-field CF 
picture yields even an exact ground state when the in- 
teraction is short-ranged enough. The cusps sticking to 
2lp,4lp, ■ ■ ■ remind us of the Tomonaga-Luttinger(TL) 
liquid, which is a totally different system in one dimen- 
sion, but the cusps, having a topological origin, do not 
change with the form of interaction. 
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FIG. 5. Full excitation spectra 
(a = 0.5, 1.0, 2.0) for TV = 9 (l F = 
ized by the gap at L — 21 f + 1(= 
SMA result is also shown (□). 



for v = 1/2 with (V 2 s-m) a 
2). The energy is normal- 
5) for each value of a. The 



The tendency that the system lies the further away 
from the Fermi liquid the shorter-ranged the interac- 
tion is consistent with analytic studies. Namely, an im- 
proved random-phase approximation (RPA) || and a 
renormalization group (RG) study ^ suggest that for 
short-ranged potential the one-particle Green's function 
has a branch cut rather than a pole just as in the TL liq- 
uid. For longer-ranged case the Fermi-liquid properties 
are recovered. To test these predictions from numerical 
low-lying spectra will require further investigations, in- 
cluding correlation function studies. However, we can 
compare the behavior of the lowest cusp (i.e., single- 
exciton branch) with the SMA, where the plm defined 
above operated on the ground state |^o) is used as the 
trial function in evaluating the energy, = {^>v\p^ LM (H — 

E )p L M\*o)/(*Q\p f LM PLM\*o) = f(L)/s{L). The SMA 
result (□ in Fig. ^|) roughly reproduces the gradient of the 
branch, although we encounter a deviation larger than 
those in the odd-fraction liquids. We can numerically 
show that the structure factor s(L) remains almost iden- 
tical as the interaction range is varied, so the change in 
the oscillator strength f(L) is dominating the shape of 
the cusp. 

To summarize, the present numerical result, done on 
the largest scale currently available, has enabled us to 
show that the gauge fluctuations in the even-fraction 
metals are substantial, but not so strong as to destroy 
the shell structure in the low-lying excitation spectrum. 
We are also extending the present study to the spin de- 
grees of freedom, which will be published elsewhere. We 
acknowledge Peter Maksym for a critical reading of the 
manuscript, and Kazuhiko Kuroki for illuminating dis- 
cussions. 
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